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IONOSPHERIC  TURBULENCE:  INTERCHANGE  INSTABILITIES 
AND  CHAOTIC  FLUID  BEHAVIOR 

I.  INTRODUCTION 

It  is  well-known  that  interchange  instabilities  produce  turbulence  in 
the  ionosphere  [Ossakow,  1979;  Fejer  and  Kelley,  1980].  The  Rayleigh- 
Taylor  instability  is  believed  to  cause  the  intense  nighttime  equatorial  F 
region  turbulence  known  as  equatorial  spread  F  (Ossakow,  1981;  Kelley  and 
McClure,  1981];  the  E  x  B  gradient  drift  instability  has  been  invoked  to 
explain  high-latitude  ionospheric  irregularities  [Keskinen  and  Ossakow, 
1983],  the  rapid  structuring  of  baritn  clouds  [Linson  and  Workman,  1970], 
and  turbulence  in  equatorial  electrojet  [Ossakow,  1979];  and  the  current 
convective  instability  may  be  responsible  for  turbulence  in  auroral 
ionosphere  [Ossakow  and  Chaturvedi,  1979a].  These  Instabilities  are 
fundamentally  similar  in  nature  in  that  they  require  a  density  gradient, 
and  they  act  to  interchange  the  high  and  low  plasma  density  regions. 
However,  they  have  different  driving  mechanisms:  the  Rayleigh-Taylor 

instability  is  driven  by  the  gravitational  force,  the  E  x  g  gradient  drift 
instability  is  driven  by  an  ambient  electric  field  or  neutral  wind,  and  the 
current  convective  instability  is  driven  by  a  current  parallel  to  the 
ambient  magnetic  field.  A  considerable  amount  of  research  has  been  devoted 
to  the  study  of  these  instabilities  and  their  application  to  ionospheric 
turbulence.  Since,  in  general,  it  is  the  nonlinear  phase  of  the 
instabilities  that  is  observed,  it  is  important  to  identify  the  salient 
characteristics  of  this  phase  (e.g.,  wave  amplitude  at  saturation,  power 
spectra,  etc.). 

Interchange  instabilities,  as  they  are  applied  to  the  ionosphere,  can 
be  divided  into  two  categories:  colllsional  and  inertial.  The  collisional 

limit  considers  »  ta  (where  v^n  is  the  ion-neutral  collision  frequency 
Manuscript  approved  September  7,  1984. 


and  «t  is  the  wave  frequency)  while  the  inertial  limit  considers  u  »  v^n» 
In  the  ionosphere,  the  transition  from  the  collisional  liait  to  the 
inertial  liait  typically  occurs  in  the  altitude  regime  *•  500  km.  The  bulk 
of  the  nonlinear  research  on  interchange  instabilities  in  the  ionosphere 
(both  theory  and  simulation)  has  been  restricted  to  the  collisional 
domain.  The  purpose  of  this  paper  is  to  extend  the  nonlinear  theory  of 
interchange  instabilities  into  the  inertial  regime. 

In  brief,  we  develop  a  set  of  mode  coupling  equations  which  describe 
the  nonlinear  evolution  of  the  Rayleigh-Taylor  and  E  x  b  gradient  drift 
instabilities.  The  model  is  restricted  to  2D  turbulence  in  the  plane 
transverse  to  the  magnetic  field,  and  only  those  modes  such  that  kL  »  1 
are  considered,  where  k  is  the  wavenumber  and  L  is  the  scale  length  of  the 
density  gradient.  He  do  not  consider  the  current  convective  instability 
for  simplicity  because  it  requires  a  component  of  the  wavenumber  parallel 
to  B  [Chaturvedi  and  Ossakow,  1979b).  He  show  that  our  mode  coupling 
equations  have  the  same  structure  as  the  equations  describing  the  Bayleigh- 
Benard  Instability;  in  particular,  when  we  consider  a  three  mode  system,  we 
show  that  the  equations  correspond  exactly  to  the  Lorenz  equations  which 
approximately  describe  the  Rayleigh-Benard  instability  [Saltzman,  1962; 
Lorenz,  1963).  Following  the  analysis  of  Lorenz  (1963),  it  is  shown  that 
the  three  mode  system  can  exhibit  a  strange  attractor  with  chaotic  behavior 
[Ruelle  and  Takens,  1971].  Ion  inertia  plays  a  critical  role  in  this 
phenomenon  in  that  if  it  is  neglected  (as  in  the  collisional  limit),  the 
three  mode  system  does  not  exhibit  chaos  and  a  stable  convection  pattern 
results. 

The  organization  of  this  letter  is  as  follows.  In  the  next  section  we 
present  the  geometry,  assumptions,  and  derivation  of  the  nonlinear  mode 
coupling  equations.  In  Section  III  we  present  our  results  for  a  three  mode 


system  and  show  the  correspondence  to  the  Rayleigh-Benard  problem.  We 
derive  a  criterion  for  the  onaet  of  chaotic  turbulence  involving  the 
dissipation  parameters  (i.e.,  diffusion  coefficient,  ion-neutral  collision 
frequency).  Finally,  in  the  last  section  we  summarize  our  findings  and 
discuss  the  implications  for  ionospheric  turbulence. 

II.  DERIVATION  OF  NONLINEAR  EQUATIONS 

The  plasma  configuration  and  assumptions  used  in  the  analysis  are 
described  as  follows.  We  consider  an  ambient  magnetic  field  in  the  z- 

A 

direction  (B  -  BQ  e^J,  a  gravitational  acceleration  in  the  -x-direction 

A 

(g  *  -g  e^),  a  density  gradient  in  the  x-direction  (n  »  and 

BUq/Bx  >0),  and  an  ambient  electric  field  in  the  y  direction 

A 

(E  -  Eg  e^J.  We  assume  two  dimensional  perturbations  in  the  xr-y  plane' 

*  A 

(transverse  to  B)  8uch  that  k  *  k  e  +  k  e  with  k_I*  »  1  and  k_L  »  1 

-  x  x  y  y  *  J 

where  L  -  [  3tn  n/3x]_l  is  the  density  gradient  scale  length.  For 

simplicity  we  assume  a  cold  ion  plasma  (T^  •  0).  We  consider  low  frequency 

turbulence  in  a  weakly  collisional  plasma  such  that  9/dt  «  0,,  vln  «  0t, 

«  0^,  and  where  ^  *  eBg/mac  is  the  cyclotron  frequency  of 

species  a,  is  the  ion-neutral  collision  frequency,  v  is  the  ion- 

electron  collision  frequency,  and  vgi  is  the  electron-ion  collision 

frequency.  We  neglect  electron-neutral  collisions  since  v  /D  «  v,  /ft, 

en  e  in  i 

in  the  F  region.  Finally,  we  consider  electrostatic  turbulence  and  assume 
quasi-neutrality  (n^  *  n^j. 

The  equations  used  in  the  analysis  are  continuity,  momentum  transfer. 


and  charge  conservation: 


0  ■  -  — -  (E  +  —  V  x  B) - -2s.-v.fv  -  V.  I 

a  v  c  -e  a_  n  el  -e  *1 

e  a 

-li.  ■  —  (E  +  —  V.  x  B)  -  v.  (V.  -  V  )  -  v.  V.  ♦  g 
at  ot  114  c  -i  ie  -i  -e  in  -i  ® 

7  •  ^  -  V  •  [nf^  -  Ye)j  -  0 


(2) 


(3) 


(4) 


We  perturb  (l)-(4>  about  an  equilibrium  and  let  n  -  n^  +  n, 

|  ■  Eg  -  7$,  and  Ya  ■  Yqq  +  Va.  To  lowest  order  in  va/na  the  equilibrium 

drifts  are  given  by  V  -  -  -  fcT  /eBl(3*n  n/3xl_l  e  (the  electron 

'*e(J  e  y 

diamagnetic  drift)  and  Vi0  -  [g/J^  +  f cEg/B)]ey  (the  ion 
gravitational  drift  and  ion  Pedersen  drift,  respectively).  Note  we  have 
chosen  VQx  ♦  VQx  -  cEQy/B.  We  solve  (2)  and  (3)  for  Ve  and  V^,  and 
substitute  these  values  into  (1)  and  (4).  We  arrive  at  the  coupled  set  of 
equations  for  a  and 


3  n 


Vn„ 


.2  n 


•  nT-  V x  •«  *  7n 


(5) 


and 


:* 

is 


v.  cE  ■»  ,  . 

in  - 1  .  a  n  _cl  (3 


57*s  +  af  r)  • 7  %•  -  f  S7  tfc  +  -r  *7i- 7  ♦  V7**  - 0 


(6) 


2 

where  De  -  i*  the 

electron 

diffusion  coefficient 

and 

we 

have 

assumed  3/3t  » 

It  may 

be  seen  that  (5) 

and 

(6) 

are 

mathematically  the  same 

as  the  Rayleigh-Benard  equations 

[Eqs. 

(17) 

and 

(18)  of  Lorenz  (1963)]  provided  the 

2 

substitution  vin  ♦  v7  i8 

made 

in  (6). 

III.  RESULTS 


K& 


We  now  present  results  of  our  analysis  for  a  three  mode  configuration. 
Prior  to  this  we  cast  (5)  and  (6)  into  dimensionless  form.  Specifically, 
we  find  that  (5)  and  (6)  can  be  written  as 


3n,  *  22  -  ~  - 

— -  ■  D  v  n  +  74 
3tx  e  1  yTl 


7n. 


7*i  x  •* 


(7) 


and 


,  3  -  -  *2  *• 

%  ♦  •«  *  7*i  •  7)7  *i 


-  v.  7  4  +  7  n. 

in  T1  y  1 


(8) 


where  tj  -  YQt,  vln  -  v^/Yq,  -  (n/n0)fL/X),  4j  -  4(c/BX2y0),  7  -  X7,  De  - 
2  1/2 

D^/X  Y0>  Yq  ■  (.(g  +  v^cEq/BJ/LJ  ,  and  X  is  half  of  the  maximum 
wavelength  permitted  (i.e.,  X  <  2L). 

We  consider  the  following  perturbations  for  4j  and  n^:  4,  *  Xj  sin  x 
sin  y  and  n^  ■  Y^  sin  x  cos  y  +  Z^  sin  2x  where  we  have  taken  1^  ■  ky  with 
kx  ■  X  *  and  ky  ■  X  1  for  simplicity.  Here,  x  and  y  represent  the  x  and  y 
spatial  coordinate  normalized  to  X.  and  only  the  coefficients  Xj.  Yp  and 
are  assumed  to  be  time  dependent.  Substituting  n^  and  4j  into  (7)  and 
(8)  we  obtain  the  following  set  of  coupled  ordinary  differential  equations: 


X  -  -  oX  +  oY 


(9) 


where  the  dot  over  a  variable  indicates  a  time  derivative,  X  *  Xj  /re,  Y  - 

Y,  /ro/2v.  ,  Z  -  rZ,,  t  -  2D  t . ,  a  -  v.  /2D  and  r  -  (4D  v.  I"1. 

1  in  l  el  in  e  e  in 

Equations  (9)-(ll)  correspond  exactly  to  the  equations  solved  by  Lorenz 

(1963)  for  the  Rayleigh-Benard  instability  (with  the  exception  that  Lorenz' 

parameter  b  is  equal  to  2  in  our  case). 

Following  Lorenz  (1963),  (9)-(ll)  can  be  analyzed  to  determine  (1)  the 

nonlinear  fixed  states  of  the  system,  and  (ii)  the  stability  of  these  fixed 

states  as  a  function  c  r  and  o.  The  fixed  states  are  determined  by  the 
•  •  • 

condition  X  -  Y  -  Z  -  0.  If  r  <  1,  the  only  stable  steady  state  is  given 

by  Xq  ■  Y0  -  Z0  -  0.  This  represents  the  state  of  no  convection,  i.e.,  the 

usual  equilibrium  state  upon  which  linear  stability  analysis  is  performed; 

in  fact,  r  -  1  is  the  point  of  marginal  linear  stability  of  the 

equilibrium.  For  r  >  1,  the  equilibritat  is  unstable  resulting  in 

convection  and  the  interchange  of  high  and  low  density  regions.  However,  r 

>  1  allows  two  additional  fixed  states,  given  by  Xq  ■  YQ  - 
1  /2 

+  [Z(r  -  1)]  and  Zq  ■  r  “  1;  these  correspond  to  convection  cells  of 

either  positive  or  negative  vorticlty.  The  linear  Instability  may  thus 
"saturate”  by  attracting  to  one  of  these  new  fixed  states  with  Xq,  Yq,  and 
Zq  providing  estimates  for  the  saturation  amplitudes  of  n  and  $.  This  is 
indeed  the  case  for  1  <  r  <  rc  =  e(o  +  5)/(o  -  3)  or  a  <  3.  In  this  range 

of  r,  the  fixed  states  are  stable  and  the  orbit  in  X,  Y,  Z  phase  space 

asymptotes  to  one  of  the  nontrivial  fixed  points.  An  example  is  shown  in 
Fig.  1  (o  ■  10,  r  ■  15)  where  the  projection  of  the  orbit  on  to  the  X-Z 
plane  is  plotted. 

For  <3  >  3  and  r  >  rc.  the  saturated  convection  pattern  described  above 
is  itself  unstable.  No  other  fixed  stable  states  exist;  this  means  that 


Che  amplitudes  X,  Y,  and  Z  oscillate  in  intensity  in  periodic  or  chaotic 
fashion.  The  magnitude  of  these  oscillations  cannot  be  determined  analyti- 
caliy  and  numerical  analysis  of  (9)  -  (11)  is  required.  Nevertheless,  some 
general  features  may  be  discerned:  even  though  the  phase  space  orbit  does 
not  approach  a  single  point,  it  does  lie  in  a  bounded,  or  "attracting", 
region  of  phase  space.  Figure  2  illustrates  such  an  orbit  for  a  -  10  and  r 
*30.  This  orbit  in  fact  tends  to  encircle  either  one  or  the  other  of  the 
two  non-trivial  fixed  states.  There  is,  however,  no  periodicity  for  the 
case  shown:  the  transition  from  encircling  one  or  the  other  fixed  points 
is  seemingly  random.  In  such  a  case  the  orbit  is  "chaotic"  and  the 

attracting  region  is  a  strange  attractor.  Since  all  orbits  are  unstable, 
one  sees  chaotic  behavior  in  the  amplitudes  of  X,  Y,  and  Z  [Ruelle  and 

Takens,  1971].  Note  that  the  amplitudes  of  X  and  Z  in  Fig.  2  can  fluctuate 
by  more  than  a  factor  of  2. 

An  important  point  to  be  recognized  concerning  the  application  of  this 
theory  to  interchange  instabilities  is  the  following.  As  noted  earlier,  we 
can  consider  two  limits:  collisional  and  inertial.  In  the  collisional 

limit  »  3/3t),  X  *  0  in  (9)  so  that  X  3  Y  and  the  problem  reduces  to 
solving  only  two  coupled  differential  equations.  For  this  case,  the 

convection  states  are  always  stable.  This  is  the  situation  that  has  been 
considered  in  most  previous  analytic  treatments  of  ionospheric  interchange 
instabilities  [Rognlien  and  Weinstock,  1974;  Chaturvedi  and  Ossakow,  1977, 
1979a, b]  although  the  stability  of  the  convection  states  was  not 
analyzed.  Thus,  for  this  simple  three  mode  system,  we  find  that  ion 
inertia  is  required  for  unstable  convection  patterns  to  occur.  We  comment 
that  the  nonlinear  behavior  of  Rayleigh-Taylor  instability  in  the  inertial 
limit  has  been  considered  by  Hudson  (1978).  However,  Hudson  (1978)  did  not 


V 


A  strange  attractpr  for  (9)-(ll)  is  illustrated  as  a 
projection  onto  the  X-Z  plane.  The  parameters  used  are  a  *  10 


?! 


J 

and  r  ■ 

30. 

Note  that  rc  - 

21.4  so  that  r  >  rc. 

All  periodic 

orbits 

are 

unstable,  as 

well  as  the  three 

equilibrium 
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find  chaotic  behavior  because  it  was  assumed  that  3/3t  ■  -  lu  where  u  is 
the  linear  eigenfrequency  in  (6).  This  ad  hoc  assumption  effectively 
reduces  the  problem  to  two  differential  equations  and  leads  to  results 
similar  to  the  collisional  limit. 

IV.  CONCLUDING  REMARKS 

We  have  shown  that  Interchange  instabilities  relevant  to  ionospheric 
turbulence  (Rayleigh-Taylor  and  E  x  B  gradient  drift),  can  be  studied  in 
the  context  of  chaotic  attractor  theory.  In  particular,  we  demonstrate 
that  for  a  simple  system  (three  modes)  the  equations  governing  these 
instabilities  are  exactly  the  same  as  those  that  govern  the  Rayleigh-Benard 
instability  [Saltzman,  1962;  Lorenz,  1963].  The  analogy  between  these  two 
instabilities  is  the  following.  The  Rayleigh-Benard  instability  is  driven 
by  a  temperature  gradient  and  convects  "hot  and  cold"  fluid  elements;  the 
temperature  gradient  is  maintained  by  a  heat  source  at  one  end  and  a  heat 
sink  at  the  other.  The  Rayleigh-Taylor  and  E  *  B  gradient  drift 
instabilities  are  driven  by  a  density  gradient  and  convect  "heavy  and 
light"  fluid  elements;  in  the  case  of  the  ionosphere,  the  density  gradient 
could  be  maintained  by  photoionization  at  one  end  and  recombination  at  the 
other.  We  have  shown  that  these  interchange  instabilities  can  exhibit  both 
stable  and  unstable  convection  patterns  in  this  system.  A  crucial  point  is 
that  unstable  convection  only  results  if  inertial  effects  are  Important, 
l.e. ,  v.  (  3/3 t. 

For  application  to  the  ionosphere,  we  consider  the  stability  of  the 
fixed  states  of  the  Rayleigh-Taylor  instability  in  the  equatorial 
ionosphere,  and  the  E  x  B  gradient  drift  instability  in  the  high  latitude 
auroral  ionosphere.  Three  observations  may  be  made.  First,  the  nonlinear 


fixed  states  given  by  Xq,  Yq,  and  Zq  correspond  to  the  saturated  potential 
and  density  fluctuation  amplitudes.  We  note  the  density  fluctuation 
amplitude  associated  with  Yq  is  n/nQ  ■  4(p#/L)(vlnv#1/YQ)1^,  while  the 
amplitude  associated  with  Zq  is  n/nQ  »  X/L.  These  estimates  agree  with 
previous  results  [Chaturvedi  and  Ossakow,  1977,  1979]  and  yield  density 
fluctuations  of  several  percent  for  typical  F  region  parameters.  These 
estimates  can  vary  by  more  than  a  factor  of  2  when  the  instabilities  are  in 
the  "chaotic  regime"  (see  Pig.  2).  Second,  since  a  •  X2vln/2velp2  u*in8 
the  normalizations  listed  after  (8),  we  find  that  103  £  a  £  10  f  for  the  F 
region  ionosphere  (200-800  km)  where  we  have  taken  X  -  500  m,  pe  ■  1.5  cm, 
vln  «  2.4  x  1011  Te1^2an  aac”1  [Strobel  and  McElroy,  1970], 
vei  "  (^ei/3«5  x  105)fne/Te3^2 )  sec”1  [Braglnskll,  1965;  Johnson,  1961] 
where  Tg  is  the  electron  temperature  in  eV,  nn  is  the  neutral  gas  density, 
ne  is  the  electron  density,  X#i  -  23.4  -  1,15  log  ne  +  3.45  log  Te,  and  the 
neutral  densities  were  obtained  from  a  Jacchla  (1975)  model  neutral  atmo¬ 
sphere.  Third,  si.nce  a  »  l,  the  critical  value  of  r  is  given  by  r£  *  o. 
For  the  Rayleigh-Taylor  instability,  the  condition  for  unstable  fixed 
states  can  thus  be  written  as  vin  <  (g/21)1^.  For  g  ■  9.8  m/sec2  and  L  - 
10  km,  we  find  that  unstable  behavior  can  occur  for  vin  <  0.02  sec-1  which 
corresponds  to  altitudes  greater  than  400-500  km  in  the  equatorial  F  region 
ionosphere.  For  the  E  x  B  gradient  drift  instability,  the  condition  for  an 
unstable  fixed  state  is  vln  <  cEq/2BL.  Taking  cEq/B  ■  6  x  102  m/sec  and  L 
■  102  km  [Weber  et  al.,  1984],  we  find  that  chaotic  behavior  can  occur 
for  vin  <  6  x  10  sec  1  which  corresponds  to  altitudes  greater  than 
roughly  500  km  in  the  high  latitude  F  region. 

We  comment  that  the  values  of  r  and  a  relevant  to  ionospheric  plasmas 
are  considerably  different  than  those  used  in  most  mathematical  studies 


(i.e.,  similar  to  those  used  in  Figs.  1  and  2).  However,  we  have  also 
performed  calculations  for  large  r  and  a  and  have  found  similar  behavior; 
namely,  the  transition  from  a  stable  attractor  to  a  strange  attractor  for 
sufficiently  large  r.  This  is  in  agreement  with  previous  studies  [Fowler 
and  McGuInness,  1982]. 

Finally,  we  have  only  considered  a  simple  three  mode  system  which 
could  be  argued  is  unrealistic.  It  is  known  that  the  nonlinear  behavior  of 
the  Lorenz  equations  for  more  than  three  modes  can  be  different  from  the 
three  mode  system  [Ott,  1981].  For  example,  we  have  shown  that  for  three 
modes,  the  nonlinear  fixed  states  are  always  stable  in  the  colllslonal 
limit  (l.e.,  viQ  >  Yq).  However,  we  have  also  developed  a  pseudo-spectral 
code  which  follows  the  evolution  of  a  many  mode  system  (~  120  modes). 
Preliminary  results  indicate  that  chaotic  fluid  behavior  can  also  occur  in 
the  colllslonal  limit,  in  sharp  contrast  to  the  simple  three  mode  result. 
A  report  on  the  nonlinear  dynamics  of  the  many  mode  system  will  follow 
shortly. 
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01 CY  ATTN  A.  JEAN  CUNCLASS  ONLY) 
01 CY  ATTN  W.  UTLAUT 
01 CY  ATTN  0.  CROMBIE 
01 CY  ATTN  L.  BERRY 

NATIONAL  OCEANIC  G  ATMOSPHERIC  AOMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULOER,  CO  80302 
01 CY  ATTN  R.  6RUBB 
01 CY  ATTN  AERONOMY  LAB  G.  REIO 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
01 CY  ATTN  I.  6ARFUNKEL 
01  CY  ATTN  T.  SALMI 
01 CY  ATTN  V.  JOSEPHSON 
01 CY  ATTN  S.  BOWER 
01 CY  ATTN  D.  OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  NA  01803 

01 CY  ATTN  RADIO  SCIENCES 

AUSTIN  RESEARCH  ASSOC.,  INC. 

1901  RUTLAND  DRIVE 
AUSTIN,  TX  78758 

01 CY  ATTN  L.  SLOAN 
01 CY  ATTN  R.  THOMPSON 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.O.  BOX  983 
8ERKELEY,  CA  94701 
01 CY  ATTN  J.  WORXMAN 
01 CY  ATTN  C.  PRETTIE 
01 CY  ATTN  S.  BRECHT 


BOEING  COMPANY,  THE 
P.O.  BOX  3707 
SEATTLE,  WA  98124 

01 CY  ATTN  6.  KEISTER 
01 CY  ATTN  D.  MURRAY 
01 CY  ATTN  6.  HALL 
01 CY  ATTN  J.  KENNEY 

CHARLES  STARK  DRAPER  LABORATORY,  INC. 
555  TECHNOLOGY  SBUARE 
CAMBRIDGE,  MA  02139 
01 CY  ATTN  D.B.  COX 
01 CY  ATTN  J.P.  6ILM0RE 

COMSAT  LABORATORIES 
LINTHICUM  ROAD 
CLARKSBURG,  MD  20734 
01 CY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 

01 CY  ATTN  D.T.  FARLEY,  JR. 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 
01 CY  ATTN  H.  L06ST0N 
01 CY  ATTN  SECURITY  (PAUL  PHILLIPS) 

EOS  TECHNOLOGIES,  INC. 

606  WilsMr*  Blvd. 

Santa  Monica,  Calif  90401 
01 CY  ATTN  C.B.  GAB8AR0 
01 CY  ATTN  R.  LELEVIER 
V.. 

ESL,  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 
01 CY  ATTN  J.  ROBERTS 
01 CY  ATTN  JAMES  MARSHALL 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KIN6  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 
01 CY  ATTN  H.H.  BORTNER 
SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.O.  80X  1122 
SYRACUSE,  NY  13201 
01 CY  ATTN  F.  REIBERT 


GENERAL  ELECTRIC  TECH  SERVICES 
CO. ,  INC. 

HNES 

COURT  STREET 
SYRACUSE,  NT  13201 

01 CY  ATTN  6.  HILLNAN 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

(ALL  CLASS  ATTN:  SECURITY  OFFICER) 
01 CY  ATTN  T.N.  DAVIS  (UNCLASS  ONLY) 
01  CY  ATTN  TECHNICAL  LIBRARY 
01 CY  ATTN  NEAL  BROWN  (UNCLASS  ONLY) 

6TE  SYLVANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  D1V 
77  A  STREET 
NEEDHAM,*  MA  02194 

01 CY  ATTN  DICK  STE1NH0F 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

01 CY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 

01 CY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA,  VA  22311 
01  CY  ATTN  J.M.  AEIN 
01 CY  ATTN  ERNEST  BAUER 
01 CY  ATTN  HANS  WOLFARD 
01  CY  ATTN  JOEL  BEN6ST0N 

INTL  TEL  8  TELE6RAPH  CORPORATION 
SOO  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

01 CY  ATTN  TECHNICAL  L18RARY 

JAYCOR 

11011  TORREYANA  ROAD 
P.O.  BOX  85154 
SAN  DIEGO,  CA  92138 

01 CY  ATTN  J.L.  SPERLING 


JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  ND  20810 

01 CY  ATTN  DOCUMENT  LIBRARIAN 
01 CY  ATTN  THOMAS  POTEMRA 
01 CY  ATTN  JOHN  OASSOULAS 

KANAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
01 CY  ATTN  T.  MEAGHER 

KANAN  TEMPO-CENTER  FOR  ADVANCED 
STUDIES 

816  STATE  STREET  (P.O  DRAWER  «0) 

SANTA  BARBARA,  CA  93102 
01CY  ATTN  DASIAC 
01 CY  ATTN  WARREN  S.  KNAPP 
01 CY  ATTN  WILLIAM  MCNAMARA 
01 CY  ATTN  B.  GAMBILL 

LINKA8IT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 
01 CY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  8  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 
01 CY  ATTN  DEPT  60-12 
01 CY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  8  SPACE  CO.,  INC. 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

01 CY  ATTN  MARTIN  WALT  DEPT  52-12 
01 CY  ATTN  W.L.  IMHOF  DEPT  52-12 
01 CY  ATTN  RICHARD  G.  JOHNSON 
DEPT  52-12 

01  CY  ATTN  J.B.  CLADIS  DEPT  52-12 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 

01 CY  ATTN  R.  HEFFNER 

M.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  MA  02173 

01 CY  ATTN  DAVID  N.  TOWLE 
01 CY  ATTN  L.  LOUGHLIN 
01 CY  ATTN  D.  CLARK 
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MCDONNEL  DOUGLAS  CORPORATION 
5301  SOLSA  AVENUE 
HUNTINGDON  BEACH,  CA  92047 
01 CY  ATTN  N.  HARRIS 
01 CY  ATTN  J.  NOULE 
01 CY  ATTN  GEORGE  NROZ 
01 CY  ATTN  W.  OLSON 
01 CY  ATTN  R.W.  HALPRIN 
01 CY  ATTN  TECHNICAL 

LIBRARY  SERVICES 

N1SSI0N  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
01 CY  ATTN  P.  FISCHER 
01 CY  ATTN  U.F.  CREVIER 
01 CY  ATTN  STEVEN  L.  GUTSCHE 
01 CY  ATTN  R.  80GUSCH 
01 CY  ATTN  R.  HENDRICK 
01 CY  ATTN  RALPH  KILB 
01 CY  ATTN  DAVE  SOULE 
01 CY  ATTN  F.  FAJEN 
01 CY  ATTN  M.  SCHEI8E 
01 CY  ATTN  CONRAD  L.  LONGHIRE 
01 CY  ATTN  8.  WHITE 
01 CY  ATTN  R.  STA6AT 

MISSION  RESEARCH  CORP. 

1720  RANDOLPH  ROAO,  S.E. 
ALBUQUERQUE,  NEU  MEXICO  87106 
01 CY  R.  STELLINGUERF 
01 CY  M.  ALME 
01 CY  L.  WRIGHT 

MITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01730 

01 CY  ATTN  JOHN  MORGANSTERN 
01 CY  ATTN  G.  HARDING 
01 CY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MAO  ISON  BLVD 
MCLEAN,  VA  22101 
01 CY  ATTN  U.  HALL 
01 CY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

01 CY  ATTN  E.C.  FIELD,  JR. 


PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

01  CY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRICS,  INC. 

4  ARROW  DRIVE 
WOBURN,  MA  01801 

01 CY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  WA  98009 

01 CY  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

R  8  D  ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY,  CA  90291 
01 CY  ATTN  FORREST  GILMORE 
01  CY  ATTN  WILLIAM  B.  WRIGHT,  JR. 

01 CY  ATTN  WILLIAM  J.  KARZAS 

01 CY  ATTN  H.  ORY 

01  CY  ATTN  C.  MACDONALD 

01 CY  ATTN  R.  TURCO 

01 CY  ATTN  L.  DeR AND 

01 CY  ATTN  W.  TSAI 

RAND  CORPORATION,  THE 
1700  MAIN  STREET  -• 

SANTA  MONICA,  CA  90406 
01 CY  ATTN  CULLEN  CRAIN 
01 CY  ATTN  ED  BEDROZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

01 CY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42nd  STREET 
NEW  YORK,  NY  10036 

01 CY  ATTN  VINCE  TRAPANI 


SCIENCE  APPLICATIONS,  INC. 

1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 

01 CY  ATTN  LEWIS  M.  LINSON 
01 CY  ATTN  DANIEL  A.  HAMLIN. 
01 CY  ATTN  E.  FRIEMAN 
01 CY  ATTN  E.A.  STRAKER 
01 CY  ATTN  CURTIS  A.  SMITH 

SCIENCE  APPLICATIONS,  INC 
1710  GOOORIDGE  DR. 

MCLEAN,  V A  22102 
01 CY  J.  COCKAYNE 
01 CY  E.  HYMAN 

SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 
01 CY  ATTN  J.  CASPER 
01 CY  ATTN  DONALD  NEILSON 
01 CY  ATTN  ALAN  BURNS 
01 CY  ATTN  6.  SMITH 
01 CY  ATTN  R.  TSUNOOA 
01 CY  ATTN  DAVID  A.  JOHNSON 
01 CY  ATTN  WALTER  G.  CHESNUT 
01 CY  ATTN  CHARLES  L.  RINO 
01 CY  ATTN  WALTER  JAYE 
01 C Y  ATTN  J.  VICKREY 
01 CY  ATTN  RAY  L.  LEAD A3R AND 
01 CY  ATTN  G.  CARPENTER 
01 CY  ATTN  G.  PRICE 
01 CY  ATTN  R.  LIVINGSTON 
01 CY  ATTN  V.  GONZALES 
01 CY  ATTN  D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

01 CY  ATTN  W.P.  BOQUIST 

TOYON  RESEARCH  CO. 

P.O.  Box  6890 
SANTA  BARBARA,  CA  93111 
01 CY  ATTN  JOHN  ISE,  JR. 

01 CY  ATTN  JOEL  6ARSARIN0 

TRW  DEFENSE  &  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90278 
01 CY  ATTN  R.  K.  PLEBUCH 
01 CY  ATTN  S.  ALTSCHULER 
01 CY  ATTN  D.  DEE 
01 CY  ATTN  0/  STOCKWELL 
SNTF/1575 


VISIDYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01803 
01 CY  ATTN  W.  REIDY 
01 CY  ATTN  J.  CARPENTER 
OlCY  ATTN  C.  HUMPHREY 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 

01 CY  ATTN:  N.  ZABUSKY 
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